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Abstract 


Detection  and  tracking  by  passive  sonar  involves  the  observation  of  the 
target-radiated  noise  in  bearing-cells  and  time-cells  that  have  their 
dimensions  predetermined  by  a  combination  of  physical  limitations,  e.g. 
the  size  of  the  acoustic  array,  and  general  factors  related  to  the  geometry 
of  the  target  motion.  This  present  study  considers  the  idealized  situation 
of  a  plane  wavefront  being  received  by  a  linear  array.  Conclusions  are 
given  relating  to  the  effects  of  target  motion  on  the  allowable  observation 
times  for  constant  speed  and  heading  conditions. 
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1.  Introduction 

The  normal  passive  detection  system  involves,  in  general,  the  assumption 
that  the  received  wavefront  of  target  radiated  noise  is  a  plane  wave.  The 
signal  processing  is  based  on  this  assumption  and  provides  a  broad-band 
'bearing  versus  time'  display  or  a  narrow-band  'frequency  versus  time' 
display.  Sometimes  the  latter  data  are  presented  as  'frequency  versus 
bearing  and  time'. 

For  a  ‘bearing  versus  time'  display  (broadband)  the  bearing  data  for  a 
given  pass-band  are  usually  obtained  from  the  outputs  of  beams  formed  from 
the  whole  array  by  integrating  them,  that  is  to  say  averaging  the  outputs 
over  some  appropriate  time  —  and  thus  the  average  represents  beam  output 
power.  The  beams  are  made  to  overlap  at  the  3-dB  down  points  so  that  the 
width  of  the  uearing  cells  can  then  taken  to  be  equal  to  the  3-dB  beamwidth 
of  the  beams. 

Alternatively,  the  bearing  data  may  be  obtained  from  the  display  of  the 
correlogram  formed  by  applying  the  split-beam  process*  to  a  given  pass- 
band.  The  cell  width  is  then  taken  as  the  effective  3-dB  width  of  this 
correlogram. 

The  time  duration  tj  of  the  cells  for  the  'bearing  versus  time'  display 
is  the  coherent  integration  time,  and  is  determined  more  by  the  cycle  time 
of  the  display  system  than  any  other  consideration. 

For  a  'frequency  versus  time'  display  {narrow  band)  he  frequency  infor¬ 
mation  in  a  given  beam  is  the  power  output  of  the  frequency  analysis  device 
applied  to  the  beam  signals.  The  effective  coherent  integration  time,  ra  , 
is  then  given  by  the  inverse  of  the  frequency  resolution,  i,e.  1/A f  (the 
frequency  cell-width).  Thus,  for  narrow-band  processing,  the  cell  size  is 
taken  to  be  Af  (in  frequency)  by  1/A f  (in  time),  where  Af  is  the  3-dB 
width  of  the  frequency  analysis  window. 

In  the  above  it  is  revealed  that  in  the  broad-band  processing  the  size  of 
the  'observation  cells'  is  fixed;  in  consequence  information  can  be  lost 
if  the  signal  moves  through  the  cells  in  a  time  less  than  the  coherent 
integration  time.  On  the  other  hand,  in  the  narrow-band  processing  some 
flexibility  is  usually  avai’able.  However  the  minimum  integration  time  is 
linked  to  the  desired  frequency  resolution,  and  so  if  the  target  signal 
moves  through  the  frequency  cell  in  a  time  less  than  this  integration  time 
then  once  again  information  can  be  lost. 

This  paper  examines  the  range  of  values  of  coherent  integration  times 
likely  to  be  available  in  ASW  encounters  based  on  the  assumption  that  a 
single  plane  wave  is  received  from  the  target,  and  that  the  receiving 
vessel  and  target  are  moving  on  fixed  courses  at  constant  speeds. 


*Horton,  C.  "Signal  processing  for  underwater  acoustic  waves",  Washington, 
G.P.O.,  McGraw  Hill,  1969. 
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2  Random  Encounter 


2.1  GEOMETRY 

The  situation  to  be  studied  is  shown  in  Fig.  1.  The  assumptions  are  that 
the  noise  radiated  by  the  target  is  received  by  the  observer  as  a  simple 
plane  wave  and  that  the  effects  of  the  motion  of  the  observer  and  target 
can  be  described  by  geometrical  changes.  This  is  effectively  converting 
the  situation  to  a  two-dimensional  far-field  problem  and  assuming  that  the 
phase  fluctuations  or  geometric  variations  of  the  wavefront,  which  occur 
due  to  its  passage  through  the  sea,  are  negligible. 


CONSTANT  SPEEDS 
CONSTANT  HEADINGS 


Wc  ,  Wr:  Relative  crossing  and  radial  speeds 
V  Target  spaed 
u  Own  ship  spaed  (observer) 
a  Target  s  aspect  angle 
ft  Target's  bearing  relative  to  observe 
6  Target’s  bow  angle 


Fig.  i.  Geometrical  relationship  of  observer  and  target . 


2.1.1  Kinematic  relationships 


From  Fig.  1  we  have  the  relative  crossing  speed  wc  and  the  radial  (or 
opening)  speed  wr  given  by 

wc  «  v  cos«  -  u  sing, 

wr  =  v  sin6  -  u  cosg. 


And  since  cos6  =  sin(g  +  6y  -  80)  and  sin5  =  cos(8  +  6j-  -  80),  where 
it/2-  (8  +  87  -  80)  is  the  angle  on  the  target's  bow  and  8  is  the  bearing 
of  the  target  relative  to  the  observer's  heading,  we  have  therefore 
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B  •  -  ~  +  eo 


T  -  e0)  -  u  sinB, 

(Eq.  1) 

T  -  80)  -  U  COSB, 

(Eq.  2) 

3  is  given  by 

(Eq.  3) 

Now  the  Doppler  shift  of  a  frequency  f  radiated  by  the  target,  as  seen  by 
the  observer,  is 


where  c  =  speed  of  sound,  i.e. 

Af  =  -  6.66  »  10"4  f  x  wr  . 


(Eq.  4) 


With  f  in  Hz,  wr  in  m/s,  and  the  speed  of  sound  taken  to  be  1500  m/s,  the 
rate  of  change  of  frequency  f  caused  by  an  outward  acceleration  along 
radial  r  is 


f  =  —  (af)  =  -  6.66  x  10‘4  f  wr 

at 


(Eq.  5) 


From  Eqs.  2  and  3  we  have 


"r  =  (jr  -  so  j  v  *in(B  +  ®t  -  e„)  -  (ST  -  80)  v  sin(e  +  eT  -  e0) 
~  ®o)u  sins  +  v  cos(B  +  By  -  80)  -  ii  cosb. 


which  gives 


•  Wr  *  •  • 

Wp  s  —  +  ®0  u  s^n6  ■  u  cos6  -  ®T  v  Sln(8  +  8y  -  80) 
r 


+  v  C0S(B  ♦  By  -  80). 


(Eq.  6) 
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2.2  STEADY-STATE  CASE 


The  steady  state  is  the  situation  usually  treated  in  target  motion  analysis 
namely  that  in  which  the  observer  and  the  target  maintain  fixed  courses  at 
constant  velocities. 


2.2.1  Steady-state  kinematics 


Under  the  steady-state  conditions,  Eqs.  3  and  7,  for  bearing  rate  and  fre¬ 
quency  rate  respectively,  reduce  to 


•  wr 
6  =  '  r1 


wc 


f  =  -  6.66  x  10-“  f  x  , 
r 


(Eq.  10) 
(Eq.  11) 


and  Eqs.  8  and  9  become 
”  2wcwr 

P  ~  -  • 

r2 

f  =  2  x  10-3  f 


x  ±L 
r2 


(Eq.  12) 
(Eq.  13) 


2.2.2  Effects  of  target  motion 


The  mathematical  expressions  for  the  integration  times  ri  and  ia  are 
required  for  the  time-domain  analysis  and  the  frequency  domain  analysis 
respectively.  They  represent  the  optimum  observation  times  for  a  moving 
target. 


Time-domain  analysis 
Time  in  beam 

Ti 


3-dB  beamwidth  , 
bearing  rate 

9 

8 


Frequency  domain  analysis 

Time  in  frequency  bln 

Ta 


bin  width 
frequency  rate 

i  fg 
f 
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To  calculate  t  we  use  the  probability  density  functions  (PDFs)  in  Eqs.  10 
and  11,  i.e.  the  respective  functions  for  g  (bearing  rate)  and  f  (rate  of 
change  of  frequency)  in  their  steady-state  conditions,  which  involves  an 
estimate  of  the  spread  of  values  of  w ,  likely  to  be  observed  in  an  ASW 
encounter.  This  is  done  by  combining  tne  PDFs  of  likely  speeds  and  direc¬ 
tions  of  target  and  tracker. 

With  the  PDFs  calculated,  and  the  frequency-cell  width  and  the  spatial 
beamwidth  of  the  towed  array  known,  the  PDFs  of  the  the  integration  times 
can  then  be  obtained. 

The  mathematics  for  the  calculation  of  the  PDF  of  T-j  in  the  time  domain 
is  developed  in  Sect.  A.l  of  Appendix  A  for  a  random  encounter,  i.e.  for  an 
encounter  in  which  there  are  no  constraints  on  the  direction  of  either  the 
target  or  the  tracker. 

The  same  procedure  is  followed  for  the  PDF  of  ra  in  the  frequency  domain. 
This  is  shown  in  Sect,  A. 2  of  Appendix  A  where  the  PDF  of  available  obser¬ 
vation  time  (xa)  is  calculated  for  the  frequency  versus  time  display. 


2.2.3  Computer  simulation  results 


From  the  technique  outlined  in  Sect.  A.l  of  Appendix  A,  the  PDF  of  the 
crossing  speed  wc  (shown  in  Figs.  2  and  3)  is  obtained  using  the  assumed 

PDFs  of  speed  for  the  ASW  frigate,  for  the  conventional  submarine,  and  for 
the  nuclear  submarine  (see  p.  5).  In  this  time  domain  analysis  we  derive 
the  PDF  of  the  bearing  rate  P  (see  Figs.  4  and  5)  and  finally  the 
integration  times  (see  Figs.  6  and  7)  for  both  conventional  and  nuclear 
submarines. 

The  PDF  of  the  available  observation  times  derived  in  the  frequency  domain 
analysis  were  calculated  with  the  method  described  in  Sect.  A. 2  of 
Appendix  A  and  the  results  are  shown  in  Figs.  8  and  9. 

Results  obtained  using  a  Montecarlo  technique  are  shown  in  Figs.  10  and  11; 
the  associated  calculation  is  presented  in  Appendix  B. 


Fig.  2 


Time  domain  analysi ■ 
crossing  speed  for  a 


£ 


Fig ,  3 .  Time  domain  analysi. 

crossing  speed  for  a 
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me  domain  analysis probability  density  function  of  bearing 
te  (B )  for  a  conventional  submarine  target. 


Time  domain  analysis:  probability  density  function  of  bearing 

rate  (B )  for  a  nuclear  submarine  target. 
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TIME  ANALYSIS 


Time  domain  analysis:  probability  density  function  of  integration 
time  (Ji)  for  a  conventional  submarine  target. 


TIME  ANALYSIS 


Time  domain  analysis:  probabil  ity  density  function  of  integration 
time  t'l)  for  a  nuclear  submarine  target. 
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FREQUENCY  RNflLYSIS 


Q 


Fig.  8.  Frequency  doma  in  analysis:  probabil  ity  density  function  of 

available  obsevation  time  (xa)  for  a  conventional  submarine 
target . 


FREQUENCY  ANRLYSIS 


SJ 


Fig.  9.  Frequency  domain  analysis:  probability  density  function  of 

available  observation  time  (TaJ  for  a  nuclear  submarine  target. 
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NONTECRRLO  RNRLYSIS  CSStO 


Fig.  10.  Montecarlo  analysis:  probability  density  function  of  the  integra¬ 
tion  time  (r )  for  a  conventional  submarine  target. 


nONTECRRLO  RNRLYSIS  CSSND 


Fig.  11.  Montecarlo  analysis:  probability  density  function  of  the  integra¬ 
tion  time  (t )  for  a  nuclear  submarine  target. 
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3.  Conclusions 

The  results  as  presented  in  Figs.  4  through  11  are  valid  for  the  case  of  a 
32-hydrophone  array  with  half-wavelength  spacing  and  a  frequency  of  300  Hz. 
They  are  considered  to  be  representative  of  a  realistic  situation,  although 
it  would  be  of  interest  to  extend  the  study  by  incorporating  a  broader 
range  of  parameters. 

From  the  PDFs  of  the  speeds  considered  in  this  study  the  results  all  show  a 
strong  dependence  on  target-observer  range,  as  one  would  expect,  since  the 
change  in  bearing  rate  of  the  target  plays  a  major  role.  At  short  ranges 
the  PDF  of  the  integration  time  is  narrow  with  the  peak  occurring  at  a 
small  value  of  time.  At  longer  ranges  the  maximum  PDF  occurs  at  a  higher 
value  of  time  and  the  spread  increases.  This  is  why  at  long  ranges  one  has 
more  time  available  for  integration,  while  at  short  ranges  the  time  for 
integration  is  limited. 

Two  methods  have  been  used  to  calculate  the  probability  distribution  of  the 
integration  time:  the  analytical  method  described  in  Appendix  A  and  the 
Montecarlo  method  described  in  Appendix  B.  The  shape  of  the  curves  in 
Figs.  10  and  11  of  the  Montecarlo  results  show  a  good  relation  to  these 
obtained  with  the  more  complicated  analytical  method  (Figs.  6  to  9). 
However,  the  analytical  method  gives,  in  addition,  mathematical  expressions 
for  the  mean  and  the  standard  deviation  of  the  available  integration  time, 
and  these  can  be  used  for  predictions  of  processing  performance. 
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(APPENDICES 

A  -  Calculation  of  the  probability  density  of 
the  available  observation  time 

B  -  Montecarlo  technique 
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Appendix  A 


CALCULATION  OF  THE  PROBABILITY  DENSITY  OF  THE  AVAILABLE  OBSERVATION  TIME 


The  general  calculation  procedure  is  to  derive  the  probability  density 
function  (PDF)  of  the  integration  time  for  the  ASW  observer  from  the 
knowledge  of  the  PDFs  of  the  observer  and  target  velocities.  In  the  random 
encounter  case  the  directions  of  movement  of  observer  and  target  are  uni¬ 
formly  distributed  in  (0,  ?n). 


A. 1  TIME  DOMAIN  ANALYSIS 

In  the  time  domain  analysis  (bearing  versus  time  display),  the  relative 
crossing  speed  for  the  random  encounter  case  is  (see  Sect.  2  of  the  main 
text ) 


wc  =  v  sina  -  u  sine 


(Eq.  Al) 


The  bearing  rate  is 
’  wc 

8  =  —  ,  (Eq.  A2) 

r 

and  the  available  integration  time  is  related  to  the  bearing  rate  by  the 
relation 


T 


(Eq.  A3) 


where  6  is  the  3-dB  beamwidth  of  the  towed -array  pattern. 


A. 1.1  Calculation  of  PDF  of  crossing  speed 


From  Eq.  Al ,  with  v  sina  =  Zj  and  u  sine  *  z2  we  have 


wc  *  *i  *  V 


(Eq.  A4) 
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Under  the  hypothesis  that  z(  and  z2  are  random  independent  variables, 
the  PDF  of  wc  from  the  theorem  of  Papoulis*  (p.  189),  is 

+  ® 

p<wc)  =  j  Pz, (wc  +  y)  Pz,(y)  dy  •  ( Eq -  A5) 

-00  *  * 

In  order  to  determine  the  PDFs  of  z{  and  z2 ,  we  define  that 
x,  =  v  and  x2  =  sina  — *■  z;  =  v  sina  , 
and  from  Papoulis*  (p.  205) 

+  +  oo 

p(zi)'  /  —  Pi  )  p2(*2  )  dXj  =  /  —  Pj  (x, )  P2  (Zt  /Xj )  dX[ , 

-  »  xi  -  “  xi 

where  p(x?)  =  p(sina)  ,  and  since  p(sina)  d(sina)  =  p(a)  da  , 

p(sina)  =  — -  =  P.L-). 

d(sina)/da  COSa 

With  a  uniformly  distributed  between  (0,  2r  ) , 

p(a)  ‘k ' 

and  therefore 

p(sina)  =  —  x  — — -  ;  [see  Papoulis*  (p.  133)]. 

2*  COSa 


From  zl  =  v  sina, 
z 

sina  =  — J-  , 
v 


and  therefore  cosa  =  (1  -  sin2a)U2  =  1  -  [1  -  (Zj/v)2]1*  . 


*Papoulis,  A.  "Probability,  random  variables,  and  stochastic  processes", 
McOraw  Hill,  New  York,  1962. 
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Hence 


P(Z)  ) 


I  p(v )  »  —  X  - i - dv  =  i 

V  /i~V~ 

oo 

—  j  . —  -  — —  dv  =  p(v  sina)  ; 

'  -  -  /v2"-  z“f 


+  00 

/  i  p(v)  x  . -  dv 

-  oo  V  /„2  .  z2 

1 

[target]. 


At  this  point  we  need  to  specify  p(v)  ,  the  analytical  form  of  the  PDF  of 
the  target  velocity. 

The  same  procedure  is  followed  for  the  observer  term  u  sing: 
x2  =  u  ,  x3  =  sing  z2  »  u  sing. 


and 


P(z2>  =  I 


^  P2(x2)  P3(x3)  dx2  =  /  ^  i-  p2(x2)  p3  (z2/x2 )  dx2. 


With  e  uniformly  distributed  between  (0,  2n  ) , 


P(*) 


and  therefore  p(sinS)  =  —  x  -JL_ 

cose 


From  zi  -  u  sing, 
z 

sine  =  ~L 
u 

and  therefore  cose  =  (1  -  sin2g)^  =  V\  -  (z2/u)2 
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V 

* 


Hence 


p(z  )  =  ;  I  p  ( u )  «  i  x  - u- - du  =  -  ;  p(u-) —  du  =  p ( u  Sine); 

.  „  U  "  /u2  .  z2  *'  -  “  /u2  -  Z2' 

[observer]. 


We  assume  a  specific  analytical  form  for  p ( u )  ,  the  PDF  of  the  observer 
(ship)  velocity. 

The  calculation  of  the  integrals  p(Zj)  and  p{z2)  for  the  cases  in  which 
the  ASW  frigate  is  the  observer  and  the  conventional  or  nuclear  submarine 
is  the  target  are  developed  in  this  report  from  given  distributions  of  the 
velocities  u,  v. 


A . 1 . 2  Integrals  calculation:  observer  (ASW  frigate) 


ASW  FRIGATE 


umj  = 

um2  =  18 


|  /  \ 

§  o  *  '  f—i  i  i  '  r  x-4  ■  ■  '  '  i — 

0  10  20  30 

SPEED  (Kn) 


um3 


20 


p(u)  is  the  PDF  of  the  observer. 


Interval :  6  <  u  <  10 


p(u)  =  mt (u-6) , 

m  =  .Q  =  0.0225  (slope). 

4 

p(u)  =  -  0.1350  +  0.225'j, 


bj  =  -  0.1350 
B,  =  0.0225. 
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Interval :  10  <  u  <  18 
p (u )  =  0.09, 


I  n terval :  18  <  _u  <  _20 

p ( u )  =  0.09  +  m3 (u-18) , 

0  -  0.09  -  0.09  _ 

m  =  - - 

20-18  2 
p(u)  =  0.9  -  0.0450u, 


0.0450, 


l  b2  =  0.09 

(  8,  =  0. 


=  0.9 

(  83  =  -  0.0450. 


The  probability  term  is  given  by 


P(z2) 


I  J  “  — PM-  du. 
"  -  - 


Since  z2  -  u  sine, 

zz  ,/u  <  1  ;  u  >  z2 . 


The  limits  of  the  integral  are  as  follows: 


1 2  =  ^  ~  J 


1  p(u)  du. 


z2  /J7T\ 


Thus  the  total  probability  is  given  by 


P(Z2)  =  Ei  ♦ 


where 


]  ,  um  1  ( b .  ♦  B  ^  ij ) 


du , 


j  u(ni  ( b  <■  Ry  j ) 

n  '  ,  "  T 

1 ,  zl 
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Final ]y ,  we  have 
for  z2  <  umi 

p(z2)  =  Fj  (z2)  +  F2(z2  =  urai )  +  F3(z2  =  um2), 

^or  umi  4  z2  <  um2 

P(z2)  =  ^2^2^  +  ^3^2  ”  uITl2  ) » 

^°r  um2  <  z2  <  Um3 

P(Z2)  =  F,(z2). 
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i 


A. 1.3  Integrals  calculation:  target  (conventional  SSK) 


> 

"a 

> 

w 

2 

UJ 

Q 


o.H 


.  0.05  i 


CQ 

< 

OQ 

O 

cc 

CL 


■  CONV.  SUBMARINE 


L_. 


i  1  - '  ■  '  \ . 

10  20 
SPEED  (kn) 


IcT 


vmi  -  3 
vm2  =  8 

vm3  =  10 
=  16 
vms  =  30 


Interval :  0  <  v  <  3 

p(v)  =  y0  +  mj (v-x0) , 

p(v)  =  y0  +  +  i^v  =  aL  +111^ 

m,  =  O-1  -  0  =  0.0333, 

1  3  -  0 

p ( v )  =  0.0333v, 


Interval :  3  <  v  <  8 


p(v)  =  a2  +  v , 
p( v )  =  0.1, 


ai  =  yo  +  m!Xo, 


/  a;  =  0 

(  Aj  =  ntj  =  0.0333. 


=  0.1 


=  m2  =  0. 
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Interval:  Aj  v_<_  10 

p ( v )  =  0.1  +  (v-H) , 


hi  . 

P  ( v ) 
P  ( v ) 


0.05  -  0.1  = 

10  -  A 

0.1  +  0.0250  v 

-0.1  +  0.0250 


=  0.0250, 

2 


-  A  x  0.0250, 


v 


I nterval  :  10  <■  v  <  In 


'  v ) 


■  \  .  r!^S  t 


a3  =  -  0.1 
A,  =  0.0250. 


^  a„  =  0.045 

(  A^  =  0. 


1 nterva 1 :  15  <  /  •  20 

p  ( v  1  =  0.045  *■  hi,  '  v  - 1 5  , 

0  -  0.045  -  0.045 

HI.  =  - - -  =  -  = 

■>  20-15  4 

p ( v )  =  0.045  -  0. 01 1 3( v-16 ) , 

p ( v )  =  0.2250  -  0.0H3v 


0.0113, 


a5  =  0.2250 
A5  =  -  0.0113. 
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total 

P(*i 

) 

is 

gi  ven 

by 

a 

v  + 

-  z2 

A 

P(Zj  ) 

=  _L 

X 

In 

mi 

_ mL 

1 

+  _L 

ir 

- 

zi 

- 

TT 

a 

V  + 

42 

-  z2 

A 

+ 

X 

In 

m2 

m2 

_ L 

+  JL 

tt 

- 

zi 

TT 

a 

v  f 

4T- 

~  I 

1  A 

+  -L 

X 

In 

m3 

vm3 

1 

+  _1 

TT 

z4 

TT 

- 

l 

_ 

a 

I\  ♦ 

45 

~2] 

A 

♦  JL 

X 

In 

mi, 

rru+ 

1 

+ 

TT 

z, 

IT 

L 

1 

_ 

a 

vm  + 

1  I 

M 

AO1 

A 

♦  -L 

IT 

X 

In 

m5 

ms 

1 

+  _L 

X 


A2  -  Z2 
mi  i 


X 


/v2  -  Z2 
m2  i 


X 


/v2  -  z2 
m3  1 


X 


42T72 

1 


X 


/v2  -  z2 
ms  1 


Since 


al 

■  0, 

h 

=  0.0333, 

vmi 

=  3, 

a2 

=  0.1, 

=  0, 

vm2 

-  8, 

a3 

it 

1 

0 

• 

►— * 

A3 

=  0.0250, 

vm3 

=  10 

a4 

=  0.045, 

\ 

■  0, 

vm4 

=  16 

a5 

=  0.2250, 

h 

=  -  0.0113, 

vms 

=  20 

and 

P(z; )  *  Fi  +  f2  *  F3  +  Fn  +  Fs  » 

we  have 

for  vmi  <  Zl  <  Vmz 

P(zt)  =  Ft  <zi )  +  f2(zi  =  %i)  +  F3(Zl  =  vmz)  +  =  vm3)  +  F5(Zj  «  vm) 
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I nterval :  6  <  v  <  16 

p ( v )  =  0.05. 

1 nterval :  16  <  v  <  30 

p ( v )  =  0.05  +  m(v-16) ,  m  =  =  -  0.3571  x  10'2. 

14 

Thus 

0  <  v  <  6,  p ( v )  =  0.8333  x  10'2  v, 

6  <  v  <  16,  p ( v )  =  0.05, 

16  <  v  <  30,  p ( v )  =  +  0.  1071  -  0.3571  x  10‘2  x  v  =  a  +  A  v, 

with 

a  =  y0  +  mx0  ;  y0  =  0.05  ,  m  =  -  0.3571  x  10-2, 

A  =  m  ;  x0  *  16  ,  a  =  0.1071. 

The  total  p(Zj)  is  given  by 
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Since 

=  0,  Ai  =  0.8333  x  10'2 ,  vmi  =  6, 

a2  =  0.05,  A2  =  0,  vm2  =  16, 

a3  =  0.1071,  A3  =  -  0.3571  x  10*2  ,  vm3  =  30, 

and 

p(Z[  )  =  F,  +  F2  +  F;j, 
we  have 
for  Zj  <  vmi 

P(2i)  ~  Fi(zi)  +  =  vmi )  +  ^  3  ( 2 1  =  vm2  )  * 

for  vmi  <  Zj  <  vm2 

Pf2l)  "  ^2  2 1  ■  *  ^3^21  ~  VIII2  )  • 
for  vm2  <  Zj  <  vm3 
P(ZX )  =  FjfZj). 


Consequently, 
for  Zj  <  vmi 
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for  vmi  <  zl  <  vm2 


a 

+  -i  x  in 


for  vffl2  <  Zj  <  vm3 


P(2 


V  + 

/  V2 

- 2 1 

m2 

m2 

1 

V 

m2 

J 

V  + 

/v2 

~2 1 

m3 

m3 

V 

I 

m2 

J 

V  + 

ft 

”1 

m3 

m3 

l 

/v2  -  V2  , 
m3  m2 


a  fv  +  -  z2l  A 

.  )  =  -1  x  in  JM - JB3 - L  +_L 

”  L  vm  J  n 
m2 


A2  -  z2  . 
m3  1 


A. 1.5  Calculation  of  PDF  of  bearing  rate 


From  the  equation 


where  wc  =  crossing  speed  and  r  =  range,  we  obtain 

•  • 

p0(e)  de  =  pWc(wc)  dwc 
or 


where 


dS  _  l 
dwc  r  ‘ 


The  PDF  of  S  is  then 
Pe(8)  =  r  pW(;(wc). 
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A. 1.6  nation  of  PDF  of  integration  time 


From  the  equation 
_  e 

e 

we  obtain 

Pt  ( f )  hi  =  p- (r  1  di 
or 

Pt  ( 1  )  -  P5  ( r  )  *  (^)_1 
(IP 


where 


dr_  =  . 

dS  (6Y 


The  PDF  of  r  is  then 


with  the  3-dB  beamwidth  given  by 

0  =  A / ( N- 1 )  d, 

in  which  a  =  wavelength, 

N  =  number  of  hydrophones, 
d  =  hydrophone  spacing. 


A. 2  FREQUENCY  DOMAIN  ANALYSIS 


In  the  frequency  domain  analysis  (frequency  versus  time  display),  the 
available  observation  time  is  given  by 


Time  in  a  frequency  bin  = 


bin  width 
frequency  rate 


(Eq.  A6) 
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Frequency  rate  is  given  by 

f  *  1-  (if)  , 

3t 


(Eq.  A7) 


where  if  is  the  frequency  Doppler  shift 


Af 


(Eq.  A8) 


in  which  wr  =  radial  velocity  component, 
c  =  speed  of  sound. 


From  Eqs.  A7  and  A8  we  have 


f  =  i.  x  f  x  wr  . 


Based  on  the  geometry  of  Sect.  2.1  of  the  main  text 


•  f  T  “c  ‘  *  * 

f  ~  c  +  9o  u  s’n6  *  u  cosS  *8IV  sin(e  +  e-p  -  e0) 

+  v  cos {0  +  eT  -  e0)J. 


(Eq.  A9) 


(Eq.  A10) 


Therefore,  from  Eqs.  Ag  and  A10  in  the  steady-state  case  we  have 


f 

c 


(Eq.  All) 


in  which  wc  =  crossing  speed  component  (m/s), 

X  =  wavelength  (m), 

r  *  range  (m). 


r  \ 
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From  Eqs.  A6  and  All  we  have 


ta  =  \  x  af0  x  r  X  1_  , 


where  Af0  is  frequency  resolution. 

Alternatively, 

fc. 


"c 


with  ka  =  *  x  Af0  x  r  . 

To  calculate  the  PDF  of  Ta,  let  xa  =  y  and 
Hence 


And  from 


we 


Py(y)  dy  *  Pz(z )  dz 


obtain  py (y)  *  Pz ( -2 )  *  I  SI'1 


d  "2ka 

But  from  Eq.  A14,  *  — —  • 

dz  z3 


(Eq.  A12) 


(Eq.  A13) 


wc  =  Z  In  Eq.  A13. 


(Eq.  A14) 


and  thus  py(y)  *  *  *3  *  Pz(z)* 

2ka 
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Now  since 


Finally,  reverting  to  our  original  notation  of  Eg.  A13,  we  obtain 


PTa(*a)  = 


.  (Eq.  A15) 


And,  if  pW(.  is  an  even  function  of  wc  , 


(Eq.  A16 ) 
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Appendix  B 


MONTECARLO  TECHNIQUE 


Tne  integration  time  t  is  a  function  of  the  following  random  variables  (see 
Sect.  ?  of  the  main  text): 


u  =  speed  of  observer, 
v  =  speed  of  target, 
u  =  bow  angle  of  target, 
p  =  bearing  relative  of  target  to  observer, 


and  is  derived  from  three  equations: 


wc  =  v  sinu  -  u  sinB,  (Eq.  Bl) 

•  w 

3  =  ~  (r  =  range),  (Eq.  B2) 

r 

t  =  —  (•;  =  3-dB  heamwidth  of  array).  (Eq.  B3) 

ti 


Therefore  the  integration  time  is  a  function  of  four  random  variables  as 
fol lows : 


i  5  F(u,  v,  a,  8 )  [r  and  e  are  constants]. 

To  calculate  the  PDF  of  t  in  terms  of  the  PDFs  of  u,  v,  a,  B  ,  the 
Montecarlo  method  is  applied  to  the  case  of  a  random  encounter  in  which  a 
and  8  are  uniformly  distributed  between  (0,  2* )  and  the  distributions  of 
u  and  v  are  given  (see  Appendix  A  for  the  assumed  distribution  of  u, 
v). 


The  Montecarlo  calculation  is  based  on  the  subroutine  RANDO  1(X),  which 
generates  random  numbers  in  the  interval  (0,  1),  and  is  available  on  the 
Uni  vac  Computer  at  SACLANTCEN. 

The  operations  on  which  the  Montecarlo  computer  simulation  is  based  are  as 
fol lows: 
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Extract  a  random  number  from 
CALL  RANDO  1{ X) 

The  output  is  a  number  between  0  to  1  and  indicates  the  interval  of  Z 
(see  Table  Bl). 

Again,  extract  a  random  number  from 
CALL  RANDO  1 ( X) 

The  output  is  again  a  number  between  0  and  1  and  will  determine  the 
value  of  the  random  variable  in  the  selected  interval  (see  Fig.  A1  in 
Appendix  A  (0  <  x  <  1)  and  Table  Bl). 

Calculate  r  correspondi ng  to  the  random  values  obtained  from 
Operation  2. 


SACLANTCEN  SM-192 


TABLE  B1 


Cumulative  interval  probability  (I)  for  observer  and  target 


(a)  ASW  frigate  (Ohserverl 


PROBABILITY 

F 

CURVE  (0  <  x  <  1) 

p(  6  <  u  •  10)  =  0. IB 

0.  IB 

u,  =  6  +  4x 

p( 10  -  u  .  IB)  =  0.72 

0.90 

u‘  =  10  +  8x 

p( 18  u  20)  =  0.10 

1 

=  18  +  2x 

(b)  Conventional  submarine  (Target) 


PROBABILITY 

z 

CURVE  (0  <  x  <  1) 

p(  0  <  v  <  3)  =  0.14 

0.14 

v  =  3x 

p(  3  <  v  <  8)  =  0.45 

0.59 

v2  =  3x 

p(  8  <  v  <  10)  =  0.12 

0.71 

v,  =  8  +  2x 

p( 10  s  v  <  14)  =  0.17 

0.88 

v,,  =  10  +  4x 

p ( 1 4  <  v  <  20)  =  0.12 

v5  =  14  +  6x 

(c)  Target  bearing  (6)  relative  to  observer 


PROBABILITY 

Z 

CURVE  (0  <  x  <  1) 

p(  0  <  B  <  2w)  *  1 

1 

S  j  =  2n  X 

(d)  Target  bow  angle  (a) 


PROBABILITY 

Z 

CURVE  (0  <  x  <  1) 

» 

p(  0  <  a  <  2w )  •  1 

1 

aj  =  2ttx 
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